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A systematic study of the impact of annealing on the electronic properties of single InAs/GaAs 
quantum dots (QDs) is presented. Single QD cathodoluminescence spectra are recorded to trace the 
evolution of one and the same QD over several steps of annealing. A substantial reduction of the 
excitonic fine-structure splitting upon annealing is observed. In addition, the binding energies of 
different excitonic complexes change dramatically. The results are compared to model calculations 
within 8-band k-p theory and the configuration interaction method, suggesting a change of electron 
and hole wave function shape and relative position. 

PACS numbers: 78.67.Hc, 73.21. La, 71.35.Pq, 78.60.Hk 



Single quantum dots (QDs) are building blocks for nu- 
merous modern devices including single-photon emitters 
[l| and storage devices 0, Q. Targeted tailoring of their 
electronic properties is of utmost importance for device 
functionality. A vanishing exciton fine-structure splitting 
(FSS) of QDs, for example, is the key parameter for gen- 
erating entangled photon pairs 0, @ . FSS post-growth 
modification was recently demonstrated using external 
electric [6j and magnetic fields [7[ as well as externally 
applied stress [|[ . In order to limit the complexity of a fi- 
nal device, it is desirable to modify the FSS permanently, 
e.g., by precise variation of the QD structure. Anneal- 
ing can considerably alter the electronic structure of QD 
ensembles leading to a change of FSS an d biex- 

citon binding energy [ll[. While Refs.Tg andlldl measured 
average properties of QD ensembles, Young et al. [ll| de- 
termined the electronic properties of individual QDs by 
performing single QD spectroscopy. However, they ana- 
lyzed QDs randomly chosen before and after annealing. 
In contrast to that, in this letter a systematic study of 
the influence of annealing on the emission characteristics 
of one and the same QD for consecutive steps of anneal- 
ing is presented. Fine-structure splittings and excitonic 
binding energies are determined. By comparing the ex- 
perimental results to calculations within the framework 
of 8-band k-p theory and the configuration interaction 
(CI) method [12| , the impact of annealing on size, posi- 
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FIG. 1: In the upper part, monochromatic CL images of the 
mesa are shown for the as-grown case and after two consec- 
utive annealing steps. The position of the examined QD is 
indicated by a white arrow. The lower part shows CL spec- 
tra of the QD before and after the two respective annealing 
steps. Lines originating from neutral and charged excitonic 
complexes are identified. 



tion and orientation of the participating wave functions 
is discussed. 

The InAs QDs were grown in an AIXTRON 200/4 
metal-organic chemical vapor deposition (MOCVD) reac- 
tor in GaAs matrix on GaAs(OOl) substrates. A 300 nm 
thick GaAs buffer layer followed by a 60 nm Alo.6Gao.4As 
diffusion barrier and 90 nm GaAs were grown. For the 
QDs nominally 1.9 monolayers of InAs were deposited 
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followed by a 540 s growth interruption. Subsequently, 
the QDs were capped with 50 nm GaAs. Finally, a 20 nm 
Alo.33Gao.67As diffusion barrier and a 10 nm GaAs cap- 
ping layer were deposited. During the growth interrup- 
tion, the QDs undergo a ripening process, i.e., material 
is transported from small to large QDs, leading to a red- 
shift of the ensemble luminescence ( 1 31 ] . Here the ensem- 
ble peak is centered at 1.06 eV at T=10 K. While most 
small QDs dissolve during the long growth interruption, 
some remain, leading to an ultralow QD density (< 10 7 
per cm 2 ) in the 1.25-1.35 eV spectral range. 

The two consecutive annealing steps lasted five min- 
utes at 710 °C and 720 °C, respectively, performed un- 
der As atmosphere in the MOCVD reactor in order to 
stabilize the sample surface. 

The sample was examined using a JEOL JSM 840 scan- 
ning electron microscope equipped with a cathodolumi- 
nescence (CL) setup. The luminescence was dispersed by 
a 0.3 m monochromator equipped with a 1200 lines/mm 
grating. The light was detected using a liquid-nitrogen 
cooled Si charge-coupled-device camera. The resolution 
of the setup is «140 [ieV. Using line shape analysis, 
the energetic position of single lines could be determined 
withtin an accuracy better than 20 /ieV. The sample was 
mounted onto the tip of a helium-flow cryostat provid- 
ing temperatures as low as 6 K used for all experiments 
throughout this letter. 

Circular mesas of 24 /im in diameter were etched into 
the sample surface. Figure [T] shows monochromatic CL 
images of such a mesa viewed from the top. One can 
clearly see, that the QD density is low enough to identify 
individual QDs and their position within the mesa. It is 
thus possible to relocate one specific QD after the sample 
has undergone an annealing step. 

The QD spectra change dramatically due to a modifi- 
cation of the QD structure during the annealing proce- 
dure. Figure shows the effect of the annealing steps on 
the spectrum of a particular QD. Neutral excitons (X), 
biexcitons (XX) and charged [positively (X + ), negatively 
(X — )] excitons could be identified following Ref. LL8l For 
the QD shown in Fig. [5J the X transition energy increased 
by 43.6 meV after both annealing steps . For ease of com- 
parison of the annealing induced variations, the energetic 
position of the X line has been set to meV. The XX 
shifts to lower energies with respect to the X line, chang- 
ing its character from anti-binding (E xx = — 2.1 meV) 
to binding (E xx = 2.6 meV) with a total change in 
binding energy AE^ X = 4.7 meV. Likewise, the X+ 
binding energy increases by AE X+ = 6.3 meV. The X~, 
on the other hand, shows the opposite trend, becom- 
ing less binding with its binding energy decreasing by 
AE X _ = -1.3 meV. The annealing process thus has a 
surprisingly large impact on all binding energies, sug- 
gesting a drastic change of the involved wave functions 
and/or their mutual position. 

Additionally, the excitonic FSS was recorded by per- 
forming polarization dependent measurements (Fig. [3]). 
For this particular QD it decreased from 170 /ieV to less 
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FIG. 2: Effect of two consecutive annealing steps on the spec- 
trum of a single QD. meV corresponds to the respective neu- 
tral exciton recombination energy (1.2738 eV for as-grown, 
1.3002 for 710 °C, and 1.3174 eV for 720 °C). The X~-line 
shifts to higher energies with respect to the X-line (i.e. the 
X - becomes less binding). All other lines follow the opposite 
trend (i.e. become more binding). 



than 20 fj,eV. It should thus be possible to reduce the FSS 
in a controlled way below the homogeneous linewidth of 
the X-line (which is on the order of a few ^eV at liq- 
uid He temperatures), a prerequisite for the generation 
of polarization-entangled photon pairs The general 

trend of decreasing the FSS [{J [l0|, [ll[ and increasing the 
XX binding energy [ll| by annealing has also been ob- 
served by other authors. However, only by measuring the 
same QD before and after annealing a detailed analysis of 
the interplay of QD morphology, wave function position 
and shape, and excitonic properties can be conducted. 
While in this letter we only show results for one specific 
QD, the described observations and trends are typical for 
all analyzed QDs. 

The experimental results were modeled using 8-band 
kp theory for the single particle orbitals and the CI 
method for the few-particle states. The as-grown QD was 
assumed to have truncated pyramidal shape with {101} 
side facets, a height of 1.42 nm, a lateral size of 11.3 nm, 
and an In content of 100 %. There is some uncertainty 
concerning these numbers since the determination of the 
structure of a few small QDs within an ensemble of large 
QDs is very difficult. While reliable information about 
the structure of similarly grown QDs exists [13j , the in- 
fluence of the long growth interruption on the morphol- 
ogy of the QDs is unknown. Our model QD yields an 
X transition energy of 1.08 eV. In our experiments it is 
likely that a slight fraction of Ga is incorporated into the 
QDs during growth, leading to a higher emission energy 
than for the model QD. The annealing process and the 
resulting exchange of In and Ga atoms grade the inter- 
faces between matrix material and QD. This effect was 
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FIG. 3: Polarization dependent measurements of the same 
QD depicted in Fig.[2]reveal the excitonic fine-structure split- 
ting. It decreases from 170 ± 20 /j,eV to less than 20 /ieV after 
the second annealing step. 

simulated by applying a smoothing algorithm for each an- 
nealing step, corresponding to Fickian diffusion. The re- 
sulting atomistic structures were relaxed with a recently 
developed bond-order potential of the Abell-Tersoff type 
that is particularly suited for InAs/GaAs heterostruc- 
tures [14j . To be more precise, we performed slab calcu- 
lations with periodic boundary conditions in the lateral 
directions, a conjugate-gradient scheme to minimize the 
total energy until the maximum force on an atom in the 
system was below 1 meV/A, and the scheme outlined in 
Ref. [HI to determine the atomistic strain tensor. The 
single-particle orbitals were then derived from a strain- 
dependent 8-band kp Hamiltonian accounting for band 
coupling and band mixing. The piezoelectric field was 
computed using the first and second order piezoelectric 
tensor [l7| . The excitonic states were calculated us- 
ing the CI method by expanding the respective excitonic 
Hamiltonian into a basis of antisymmetrized products of 
single particle wave functions built from a total of six 
electron and six hole wave functions. This procedure ac- 
counts for direct Coulomb interaction, exchange, and in 
part correlation [l8l |. 

Results for the binding energies of X + , X~ , and XX 
are shown in Fig. 2] The experimentally observed en- 
ergetic order of X~ and X + is well reproduced by our 
theory. Also, the general trend of the binding energies 
upon annealing (X~ becomes less binding, X + becomes 
more binding) is in good agreement. However, there 
are some discrepancies between theory and experiment: 
The change of character of the biexciton (i.e. from anti- 
binding to binding) as observed in the experiment is not 
reproduced by the modeling. Also, the experiment re- 
vealed a much stronger relative shift for X + and XX than 
for X- (see Fig.©. 

These discrepancies imply that an important effect is 
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FIG. 4: Calculated binding energies for X - , X + , and XX. 
Open symbols correspond to a reduction of electron-hole 
Coulomb interaction by 0.6 meV simulating a misalignment 
of electron and hole wave functions (see text). 

not correctly modeled. Remarkably the disagreements 
can be resolved, when the attractive Coulomb interaction 
between electron and hole (J e h) is artificially lowered for 
the as-grown QD by 0.6 meV, corresponding to a reduc- 
tion of 3%, while repulsive electron-electron (J ee ) and 
hole-hole (Jhh) Coulomb terms are left unaltered. Such 
an effect may result from a reduced electron-hole wave 
function overlap caused by either a slight elongation and 
misalignment of the wave functions or a mutual vertical 
displacement. 

Using the assumption of reduced J e h for the as-grown 
QD we observe both the biexciton crossover as well as the 
relative insensitivity of X~ binding energy to annealing 
as seen in the experiment (open symbols in Fig. 0]). A 
symmetrization of the wave functions with annealing is 
further supported by the drastic reduction of FSS [§]. 

In conclusion, we have recorded emission spectra of 
single QDs and followed their evolution under an anneal- 
ing procedure. The binding energies of different excitonic 
complexes change on the order of several meV and the 
FSS decreases from 170 /ueV to less than 20 fieV. These 
results can be understood by a change of electron and 
hole wave function shape and mutual position. We have 
thus demonstrated a powerful tool to tailor single QDs' 
electronic properties for their use in potential applica- 
tions. In particular, zero FSS, essential for emission of 
entangled photon pairs, can be enforced. 

The calculations were performed on the IBM 
pSeries 690 supercomputer at HLRN within project 
No. bep00014. This work was supported by the Deutsche 
Forschungsgemeinschft in the framework of SfB 296 and 
the SANDiE Network of Excellence of the European 
Commision, Contract No. NMP4-CT-2004-500101. 



4 



[1] O. Benson, C. Santori, M. Pelton, and Y. Yamamoto, 

Phys. Rev. Lett. 84, 2513 (2000). 
[2] S. Cortcz, O. Krebs, S. Laurent, M. Senes, X. Marie, P. 

Voisin, R. Ferreira, G. Bastard, J.-M. Gerard, and T. 

Amand, Phys. Rev. Lett. 89, 207401 (2002). 
[3] A. Marent, M. Geller, D. Bimberg, A. P. Vasi'ev, E. S. 

Semenova, A. E. Zhukov, and V. M. Ustinov, Appl. Phys. 

Lett. 89, 072103 (2006). 
[4] C. Santori, D. Fattal, M. Pelton, G. S. Solomon, and Y. 

Yamamoto, Phys. Rev. B 66, 045308 (2002). 
[5] R. M. Stevenson, R. M. Thompson, A. J. Shields, I. Far- 

rer, B. E. Kardynal, D. A. Ritchie, and M. Pepper, Phys. 

Rev. B 66, 081302(R) (2002). 
[6] K. Kowalik, O. Krebs, A. Lemaitre, S. Laurent, P. Senel- 

lart, P. Voisin, and J. A. Gaj, Appl. Phys. Lett. 86, 

041907 (2005). 

[7] R. M. Stevenson, R. J. Young, P. See, D. G. Gevaux, 

K. Cooper, P. Atkinson, I. Farrer, D. A. Ritchie, and A. 

Shields, Phys. Rev. B 73, 033306 (2006). 
[8] S. Seidl, M. Kroner, A. Hogele, K. Karrai, R. J. Warbur- 

ton, A. Badolato, and P. M. Petroff, Appl. Phys. Lett. 

88, 203113 (2006). 
[9] W. Langbein, P. Borri, U. Woggon, V. Stavarache, D. 

Reuter, and A. D. Wieck, Phys. Rev. B 69, 161301(R) 

(2004). 

[10] A. I. Tartakovskii, M. N. Makhonin, I. R. Sellers, J. 



Cahill, A. D. Andreev, D. M. Whittaker, J.-P. R. Wells, 
A. M. Fox, D. J. Mowbray, M. S. Skolnick, K. M. groom, 
M. J. Steer, H. Y. Liu, and M. Hopkinson, Phys. Rev. B 
70, 193303 (2004). 
[11] R. J. Young, R. M. Stevenson, A. J. Shields, P. Atkin- 
son, K. Cooper, D. A. Ritchie, K. M. Groom, A. I. Tar- 
takovskii, and M. S. Skolnick, Phys. Rev. B 72, 113305 
(2005). 

[12] O. Stier, M. Grundmann, and D. Bimberg, Phys. Rev. B 

59, 5688 (1999). 
[13] U. W. Pohl, K. Potschke, A. Schliwa, F. Guffarth, D. 

Bimberg, N. D. Zakharov, P. Werner, M. B. Lifshits, 

V. A. Shchukin, and D. E. Jesson, Phys. Rev. B 72, 

245332 (2005). 

[14] T. Hammerschmidt, P. Kratzer, and M. Scheffler, unpub- 
lished. 

[15] C. Pryor, J. Kim, L. W. Wang, A. J. Williamson, and A. 

Zunger, J. Appl. Phys. 83, 2548 (1998). 
[16] G. Bester, X. Wu, D. Vanderbilt, and A. Zunger, Phys. 

Rev. Lett. 96, 187602 (2006). 
[17] G. Bester, A. Zunger, X. Wu, and D. Vanderbilt, Phys. 

Rev. B 74, 081305 (2006). 
[18] S. Rodt, A. Schliwa, K. Potschke, F. Guffarth, and D. 

Bimberg, Phys. Rev. B 71, 155325 (2005). 



